Manipulation of Giant Faraday Rotation in Graphene Metasurfaces 
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Faraday rotation is a fundamental magneto-optical phenomenon used in various optical control and magnetic 
field sensing techniques. Recently, it was shown that a giant Faraday rotation can be achieved in the low-THz 
regime by a single monoatomic graphene layer. Here, we demonstrate that this exceptional property can be 
manipulated through adequate nano-patterning, notably achieving giant rotation up to 6THz with features no 
smaller than lOOnm. The effect of the periodic patterning on the Faraday rotation is predicted by a simple 
physical model, which is then verified and refined through accurate full-wave simulations. 



Thanks mainly to its unconventional band structure and 
two-dimensional nature |fl] 13, graphene exhibits many ex- 
ceptional electrical and optical properties such as high elec- 
tron mobility (3]|4), quantum Hall effect [5], highly confined 
plasmonic propagation [6 1, large nonlinear Kerr effect Q and 
giant Faraday rotation QQ. Moreover, these properties can 
be dynamically tuned through electrostatic and magnetostatic 
bias, which is another major feature of graphene. In particular, 
Faraday rotation, which is the focus of this work, is used in 
various applications in optical telecommunications and laser 
technology (modulators, optical isolators and circulators) J5J, 
as well as for advanced magnetic field sensing and even in 
spintronics [ 10 1. 

Nevertheless, despite the tuning possibility, the general be- 
havior of graphene devices is determined by the intrinsic char- 
acteristics of the material and can hardly be tailored for a spe- 
cific application or wavelength of interest. In the case of the 
giant Faraday rotation, a 7T magnetic bias allows the dynamic 
differential rotation of a plane wave polarization up to about 
6 degrees. However, the rotation is maximum only at low fre- 
quencies and decreases for smaller wavelengths, notably re- 
ducing very rapidly above 2THz fH. In ifTTIl . a scheme based 
on electric biasing of graphene and consequently increasing 
the chemical potential is proposed to tailor the Faraday rota- 
tion property of graphene for microwave applications. 

Following the advent of the metamaterial concept, it was 
shown that patterning bulk materials to control their electro- 
magnetic properties enables many fantastic applications such 
as negative refraction, photonic bandgap and perfect lenses 
|[T2l IT3l . This idea was then applied to graphene in lfl4l for 
the control of the transmission magnitude, and other func- 
tionalities such as absorbers and cloaks were also recently 
suggested in |[T5l . Micro-patterning of graphene monolayers 
can be achieved for instance using electron beam lithography 
(EBL) and oxygen plasma etching [14|, which confirms the 
applicability of graphene micro-patterning for optical devices. 
Here, we theoretically demonstrate that one of the most amaz- 
ing properties of graphene, namely, the single-atomic layer 
giant Faraday rotation, can be efficiently manipulated through 
nano-patterning. It is notably predicted that giant Faraday ro- 
tation can be achieved at much higher frequencies than us- 
ing a uniform unpatterned layer. Our approach consists in 
first comparing our numerical results with formerly measured 



results for uniform layers, and then predicting the effect of 
nano-patteming on the Faraday response. The results are also 
shown to be in good agreement with a proposed simple phys- 
ical circuit model. 

Consider first a uniform graphene layer on an infinitely 
thick SiC substrate (Fig.[T]) with Dirac voltage V d = 0.34eV, 
temperature T = 5K and phenomenological scattering rate 
r = lOmeV (8|. The frequency-dependent Faraday rotation 
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FIG. 1: Illustration of Faraday rotation and the unit cells of both 
uniform and patterned graphene layers 

numerically computed using Kubo-formalism fl6l . as well as 
the corresponding measurements (8), are shown in Fig. [2] for 
different magnetic bias fields. The results are in fair agree- 
ment, the discrepancies being due to defects in the fabricated 
graphene layer IfTTIl . 

The Faraday rotation can be well approximated as 



9(w,fl)~Zo/(M)9l{a (<M)}, 



(1) 



where Zq is the impedance of vacuum, /(co) depends on the 
substrate properties and o„(co,Z?) is the off-diagonal term in 
the conductivity tensor of graphene (a = [a c / g ;—g a^]), 
itself calculated using Kubo's formalism Graphene 
anisotropic conductivity closely follows a Drude conductivity 
model at THz. Hence, both diagonal and off-diagonal terms 
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FIG. 2: Rotation angle of the transmitted plane wave versus fre- 
quency for a uniform graphene monolayer residing on SiC substrate 
under different magnetostatic biases, experimental (see Ref. Q3) and 
theoretical 



of the corresponding surface impedance can be accurately ex- 
pressed as a superposition of resistive and inductive compo- 
nents. Indeed, by computing the impedance matrix 



Z = a- 1 = 



Z d 
Z n 



-Z 



(2) 



one obtains the impedance values shown in Fig. [3] The plots 
reveal that both the diagonal and off-diagonal terms are very- 
well approximated by an equivalent surface R-L impedance 
model, i.e. Z ( / = Rd + JCOLj and Z () = R + j(£>L () with Rd ~ 
380Q/D, L d ~ 25.5 pH/D, R a ~ —87.1 • BQ./D, and L ~ 
2.05 x 10~ 3 BpH/D. It is observed that the off-diagonal 
terms depend on the magnitude of the bias magnetic field B 
(in Tesla), which is obviously the origin of the possibility for 
dynamically controlling Faraday Rotation through magnetic 
biasing. 

Following this circuit model representation, we inferred 
that artificially creating some electrical energy reactive field 
storage along the graphene layer would allow further control 
of o„((B,B). From a practical perspective and inspired by the 
metamaterial concept, the simplest way to achieve the desired 
capacitive effect is to periodically pattern out narrow strips of 
graphene. Indeed, consider now graphene periodically pat- 
terned as in Fig.[T] The stored energy associated with the 
electric field in the gaps translates in an additional isotropic 
capacitive effect. In addition, the overall relative surface of 
graphene is reduced and thus the inherent resistive and induc- 
tive graphene impedances scale accordingly. Therefore, the 
total effective impedance can be written as follows: 
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Now, the effective capacitance associated with the periodic 
pattern considered here can be estimated as |[T8l 



e,. + l nd 

C = e L ln(csc — ), 

7C 2L 



(4) 
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FIG. 3: (a) real part of Z,/, (b) imaginary part of Zj, (c) real part of 
Z„ and (d) imaginary part of Z fl in the anisotropic impedance matrix 
of a magnetically biased graphene monolayer 



where e r is the dielectric permittivity of the substrate. Al- 
though this equation was derived for perfect conductors, the 
high electron mobility in graphene justifies this assumption 
for evaluating the distributed capacitive effect here. 

The proposed model was used to compute the effective off- 
diagonal conductivity of a patterned graphene monolayer with 
different lattice constants and B = 7T. The same parame- 
ter was also computed using an in-house full-wave simulator 
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based on the periodic Method of Moment |[T9l . which is able 
to handle graphene anisotropic conductivity as well as Flo- 
quets periodic boundaries to rigorously account for the cou- 
pling between unit cells. Note that the results of the code 
exactly fit the analytical formula available in the case of uni- 
form graphene layer (Fig.|2|. Fig.|4] shows the comparison 
of both modeling methods for two representative cases. The 
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FIG. 4: Real part of the effective off-diagonal conductivity, proposed 
physical model versus accurate full-wave approach 

good agreement demonstrates the validity of the above physi- 
cal interpretation, as well as the ability of the simple proposed 
model to provide good prediction of the frequency and bias 
dependent Faraday rotation. 

The corresponding Faraday rotations are shown in Fig. [5] 
The graphs again include the above physically-based approx- 
imate model compared with the detailed full-wave solution. 
These results clearly validate the physical interpretation and 
model provided, as well as the possibility for manipulation 
of the Faraday rotation in terms of frequency and magnitude. 
The highest frequency of giant Faraday rotation in Fig. [5] is 
between 4THz and 6THz, a region where that of the uniform 
graphene is much weaker and strongly frequency-dependent. 
In fact, this case corresponds to a gap width d = lOOnm that 
is easily achievable via e-beam lithography, thereby demon- 
strating the practical feasibility of the concept to such (and 
even smaller) wavelengths. Comparison between Fig. [5] also 
reveals that the amount of maximum giant Faraday rotation 
mostly depends on the overall relative graphene area in the 
patterned surface (this value is constant in Fig. [5^), while the 
frequency at which the maximum rotation is obtained is con- 
trolled by the absolute value of L and d. 

The reported phenomenon provides a prominent degree of 
freedom for applications of giant Faraday rotation and the as- 
sociated magneto-optical Kerr effects. Through the patterning 
of a graphene layer, a dynamically-controllable giant Faraday 
rotation can be obtained at a frequency determined a priori 
by the patterning dimensions. This technique can potentially 



lead to a variety of applications such as ultrathin circulators 
and isolators for both infrared and THz, miniaturized Fara- 
day rotators for amplitude modulation, or integrated magnetic 
field sensors. From a more general perspective, these results 
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FIG. 5: (a) Rotation angle of patterned graphene for different lat- 
tice constants and 7T magnetic bias, with gap length fixed to one- 
eleventh of the lattice constant, i.e. d = L/l 1 and (b) Rotation angle 
of patterned graphene for different gap widths and 7 T magnetic bias 
with the lattice constant fixed to L = 1 . lf/m 



motivate the further study of graphene patterning for a much 
better and versatile control of numerous exceptional electro- 
magnetic properties. 
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